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o : 

We calculate the correlation energy of a two-dimensional homogeneous electron gas using several 
available approximations for the exchange-correlation kernel / xc (g,w) entering the linear dielectric 
response of the system. As in the previous work of Lein et al. [Phys. Rev. B 67, 13431 (2000)] 
on the three-dimensional electron gas, we give attention to the relative roles of the wave number 
and frequency dependence of the kernel and analyze the correlation energy in terms of contributions 
from the (q, iui) plane. We find that consistency of the kernel with the electron-pair distribution 
function is important and in this case the nonlocality of the kernel in time is of minor importance, 
as far as the correlation energy is concerned. We also show that, and explain why, the popular 
Adiabatic Local Density Approximation performs much better in the two-dimensional case than in 
the three-dimensional one. 
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I. INTRODUCTION 



Fluids of electronic carriers with essentially two-dimensional (2D) dynamics in semiconductor structures present a 
rich phenomenology, in which dynamical correlations take increasing importance as the electron density is lowered 1 . 
Many of the electron-electron interaction effects can be understood with the help of the homogeneous electron gas 
(EG) models. A central role in the theory of short-range correlations in the EG is played by the wave number and 
frequency dependent exchange-correlation kernels or equivalently by the local field factors entering the linear response 
properties of the modeA These provide key inputs for some applications of density functional theory* and for studies 
of quasi-particle properties such as the effective mass and the effective Lande g-factor— . These properties are known 
I ' from experiment for carriers in semiconductor structures over a wide range of carrier density^. A great deal of accurate 
, information has come from quantum Monte Carlo (QMC) studies^, but theoretical understanding continues to draw 
\Q 1 interest. 

The exchange-correlation kernel / xc (9, w) entering the dielectric response of the EG determines its correlation energy 
through the fluctuation-dissipation theorem relating the imaginary part of the inverse dielectric function to van Hove's 
dynamic structure factor S(q, oS). Starting from S(q, u>) an integration over frequency is needed to obtain the electron- 
pair structure factor S(q), which in Fourier transform describes the shape of the instantaneous exchange-correlation 
(Pauli-Coulomb) hole surrounding an average electron located at the origin. From S(q) the ground-state energy E g 
can be calculated by means of further integrations over wave number and over coupling strength (see for instance Ref. 
0) and the correlation energy e c is defined, as usual, as the difference between E g and its Hartree-Fock value. 

For the 3D EG Lein et al& have analyzed in detail this procedure for calculating the correlation energy, with 
I ■ particular attention to the roles of the wave number and frequency dependence of / xc (g,w). While the wave number 
dependence of the kernel is physically related to the spatial shape of the exchange-correlation hole, as already remarked, 
its frequency dependence reflects the inertia of the hole as the electron at the origin moves through the EG. When 
the frequency dependence of f xc (q,uj) is omitted and its wave number dependence is approximated by its leading 
long-wavelength term, one is accounting for short-range exchange and correlation effects in the EG only through the 
compressibility sum rule (see for instance Ref. 0). This approximation provides the basis for the so-called Adiabatic 
Local Density Approximation (ALDA), which has often been used in dealing with time-dependent phenomena in 
jj] ■ inhomogeneous electronic systems in the low-frequency regime^. 

In the present work we extend the analysis of Lein et al& to the 2D EG by examining how several known forms of 
f xc (q,u>) perform in regard to the calculation of its correlation energy over a range of values for the coupling strength. 
The standard of comparison is the parametrized form reported for the correlation energy by Rapisarda and Senatore 9 
from their Diffusion QMC results. The paper is organized as follows. In Sec. II we report the detailed expressions 
relating e c to / xc (g,w) and in Sec. Ill we list the various forms of f xc (q,u) that we have considered. Section IV 
presents our main numerical results and discusses them with special emphasis on the usefulness of the ALDA in 2D as 
compared with 3D. Finally, in an Appendix we briefly comment on the scaling properties of the exchange-correlation 
kernel for inhomogeneous electron fluids of arbitrary dimensionality. 
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II. CORRELATION ENERGY FROM EXCHANGE-CORRELATION KERNELS 

We consider a fluid of electrons moving in a plane and interacting by the e 2 /r law. The correlation energy (per 
electron) can be written in terms of the difference between the potential energy of the fluid and the exchange energy, 
integrated over the coupling-strength parameter A: 



dX 



2 Jo A J (2tt) 



d 2 q 



^2 q 



v X q [S x {q)-SMq)]. 



(1) 



Here u A = 2nX/q is the 2D Fourier transform of the Coulomb potential with strength A, S\(q) is the structure factor 
of an EG with interactions and Shf(<?) is the Hartree-Fock structure factor, 



ShfO?) = 



arcsin (q) + q y/l — q 2 for q < 1 

for q > 1 



(2) 



where q — q/(2kp) with kp = (27m) 1 / 2 being the Fermi wave number determined by the 2D electron density n. 
The fluctuation-dissipation theorem relates the structure factor to the density-density response function (g, iit) 



calculated on the imaginary frequency axis, 



Sx(q) = -— / 
Jo 



(3) 



The response function can in turn be expressed through the frequency-dependent exchange-correlation kernel / xc (<7j w ) 
on the imaginary frequency axis, 



1 - K + f2c(q,iu)]xo(q,iu) 



where 



Xo(q, iu ) 



m 

Tift 2 
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l j-_\Jf(q,u) + ^/f 2 (q,u)+4q 2 ^ 



(4) 



(5) 



is the density-density response function of the 2D ideal Fermi gaa^fi. Here, u = muj (Hkp q) and f(q, u) — q 2 — u 2 — 1. 



From Eqs. (3) and (4) it is possible to rewrite Eq. (1) in the form 
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qdq 
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X q 1 - [v% + /4(g, iu)]xo(q, iu) 



(6) 



In the paramagnetic case of present interest e c depends only on n, or equivalently on the dimensionless density 
parameter r s defined by r s as — (7rn) -1 / 2 , as being the Bohr radius. 

In fact, the full dependence of f£ c (q,iu) on A is not needed for calculating e c at any given r s . It suffices to know 
the r s dependence of the exchange-correlation kernel at full coupling strength, i.e. f xc (q,iu) = /x C =e (q,iu)- That is, 
the correlation energy can be calculated from 



e c (r s ) 



4 I 'M)<K 
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Ryd 



where 



lc(r s ) = ^- r dq[S(q)-S HF (q)} = f°° dq H 'du^f 

kf Jo Tmk F J J Q 1 - [v. 



[Xo{q, iu)] 2 [v q + f xc (q, iu)} 



- q + fxc(q,iu)]xo{q,iu) ' 

Several known approximations to f xc (q, iu) as needed for this calculation will be recalled in the next Section. 
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III. APPROXIMATE KERNELS FOR THE HOMOGENEOUS ELECTRON GAS 



We have tested five main forms of the kernel f xc (q,iu). In three of these (designated by the acronyms ALDA, 
DPGT, and STLS) the frequency dependence of the kernel is omitted. A dynamical kernel has been used in the 
approximations designated by QSTLS and AKA. Some details are as follows. 

(i) ALDA: for the homogeneous electron gas the adiabatic local density approximation omits both the wave number 
and the frequency dependence of the kernel. It amounts to taking into account the compressibility sum rule by 
setting 

^•■••J-ikt 1 -?) (9) 

where k is the compressibility of the EG and k that of the ideal Fermi gas, given in 2D by Kq — 7rr*/2 in units of 
a^/Ryd. The ratio kq/k is related to the correlation energy (in Ryd units) by 



K TT 8 



1 de c 



dri r s dr s 



(10) 



We evaluate this input from the parametrized form of e c (r s ) of Rapisarda and Senators^. 

(ii) DPGT: an analytical expression for the static kernel / X c(<7, 0) as a function of r s in the range < r s < 10 has 
been obtained by Davoudi et alii by fitting the available Diffusion Monte Carlo data for the local field factor in the 
2D EG™. Their expression embodies the compressibility sum rule as well as the asymptotic high-g behavior of the 
static kernel. By setting 

/ x D c pGT ( 9 ^ U )=/ x r c ( 9 ,o) (ii) 

we are evidently omitting the frequency dependence of the kernel but taking into account its " exact" dependence on 
wave number at low frequency. 

(iii) STLS: the Singwi-Tosi-Land-Sjolander approximation^ omits the frequency dependence of the kernel by re- 
lating it self-consistently to the structure factor S(q) through 



f^( q , iu ) = - / „ fc [S(| q -k|)-l]. (12) 



; . S1LS; ( ;m _ 1 /^k_q_k 

n J (2tt) 2 

This expression was justified by an analysis of kinetic equations in the presence of a time- dependent weak external 
potential and has been used rather widely in the literature. 

(iv) QSTLS: a frequency dependence of the local field factor in an STLS-type self-consistent theory has been 
included by later authors. We adopt in particular the expression^* 

fT TLS (q, iu)= l -l^ V k [5(|q - k|) - 1] (13) 



nj (2tt) 2 xo(q,iu) 

where the function xo(q> k; iu) is given bjii* 



m 



Here c = q • k/g and 



Xo(q,k;iu) = _ c (l-r). (14) 



2 y 2 + z2 _ c 2 [(y 2 +z 2_ c2)2+4cV] l/2 

r= 2? + 2^ (15) 



with y = Ikpjq and z = 2mu/hq 2 . 

(v) AKA: Atwal et aim* have recently proposed an analytic expression for the dynamical local field factor which is 
based on many-body perturbation theory and on known sum rules, in the same spirit as in the work of Richardson 
and Ashcroftii on the 3D EG. With the notations of Ref. HI 

f^ A {q,iu) = -v q [G s (q,iu) + G n (q,iu)} (16) 

where 



(17) 
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and 

n ( \ a n (u)q + b n (u)q 7 

G n (q, iu) = — — -77T7 ■ (18) 
1 + c n (u)q + d n (u)q b 

Here a S)Jl (u), 6 s ,n(w), c S)T j(u) and d s , n (u) are functions of frequency and depend on two input parameters, which 
are the contact value g(0) of the pair distribution function and the correlation energy. In our calculations we have 
replaced Eqs. (41) and (52) of Ref. for d s , n (u), which contain typographical errors, by the following expressions 
communicated to us by Dr. G. S. Atwal: 

C A \°°uj 
4M = 6[C S A0 - 1 + 5 (0)](1 + + [1 - .9(0)](1 + u*) ' (19) 



" l ; (5C„A° + 2A«)(1 + uj a ) [57„wVa + 2A~(1 + 0.38w) 1 /2](i + 0.38^(1 + w 4 ) " 1 J 

For expressions of the other functions and of the various parameters the reader is referred to the original paper. 

In determining the input for this form of the kernel we have again used the form of Rapisarda and Senatore^ for 
e c (r s ) and followed Atwal et ali& in taking the values of g(0) from an interpolation formula proposed by Polini et al~*. 
We have also checked simple variants of this scheme, taking g(0) from calculations in the ladder approximation's or 
dropping the frequency dependence of /^. KA (q, iu). 



IV. NUMERICAL RESULTS AND DISCUSSION 



We turn to a presentation of our numerical results, which are collected in Figures 1-4. In Figure 1 (top panel) 
we show the difference between the correlation energy calculated by means of Eqs. J7J and JSJ and the "exact" 
correlation energy eJ? MC from the parametrization of Ref. 0, for electron densities in the range 0.5 < r s < 10. In 2D 
this range extends into the intermediate-to-strong coupling regime. In the bottom panel of Figure 1 we report the 
calculated correlation energies in the weak-coupling regime < r s < 0.8, to show that none of the theories that we 
have examined satisfies the exactly known limiting expression^ for r s — > 0, 

2 \X2 

e c (r s ) -» -0.385 - -^-(10 - 3vr)r s lnr s + o(r s ) . (21) 
Sir 

This panel also reports the correlation energy given by the Random Phase Approximation (RPA), which corresponds 
to setting the exchange-correlation kernel to zero and evidently is badly in error even at low r^. 

Turning to the results shown in the top panel of Figure 1, a pleasant surprise in 2D is that ALDA does not 
overestimate badly the correlation energy as was shown by Lein et al& to be the case in 3D. In fact, in the range 
2 < r s < 10 the ALDA results are seen to be quite close to the DPGT results, in which the full nonlocality of the 
static kernel in space has been included. The reason for this can be inferred from the work of Davoudi et alml, 
who showed that in 2D satisfying the compressibility sum rule reproduces quite accurately the g-dependence of the 
static kernel almost up to q ~ 2kp. Following Lein et al& we have also examined frequency-dependent forms of the 
long- wavelength kernel and in particular that proposed by Qian and Vignale 2 — . We have found that this overestimates 
the correlation energy by as much as ~ 0.1 Ryd at r s = 3. It appears, therefore, that nonlocality of the kernel in 
space should be included if its nonlocality in time is accounted for in the calculation of the correlation energy. 

This is, of course, what is being done in the AKA calculation. The average deviation of its results in the top panel 
of Figure 1 is approximately 0.02 Ryd, which is somewhat better than in ALDA or in DPGT. Taking g(0) in the 
AKA kernel from the ladder-diagram calculation^ raises this deviation to about 0.03 Ryd and, most importantly, 
dropping the frequency dependence of the AKA kernel brings it to about 0.04 Ryd. The summary conclusion of all 
these calculations is, therefore, that as far as the correlation energy in 2D is concerned one does already fairly well 
by taking a constant exchange-correlation kernel adjusted to the compressibility sum rule and that nonlocality of the 
kernel in both space and time needs including in order to obtain some improvement in the results. 

Figure 1 in its top panel further shows that STLS gives a very good estimate of the correlation energy, with a 
deviation of about 0.004 Ryd on average. In this case inclusion of the frequency dependence of the kernel by means of 
the QSTLS recipe has only a minor effect, bringing the average deviation down to about 0.0035. Consistency of the 
exchange-correlation kernel with the electron-pair structure thus appears to be important, as far as the calculation 
of the correlation energy is concerned. This fact may give a useful suggestion for further improvement of the more 



5 



sophisticated kernels that have been proposed in the more recent work. We may also remark that Dobson et aim* 
have recently extended the STLS scheme to inhomogcncous electronic systems and shown that it gives good results 
for the correlation energy of jellium slabs of finite thickness, lying within 3% of the Diffusion Monte Carlo results of 
Acioli and CeperleyS^. 

We turn now to Figures 2-4, in which following again Lein et al& we provide an analysis of the correlation energy 
at r s = 1 into contributions from correlations between density fluctuations of different wave vectors and different 
imaginary frequencies. Equations Q and (JSJ) naturally define a wave- vector analysis e c {q) if we write 

/>oc 

£c= / e c (q)d(q/k F ) (22) 
Jo 

and an imaginary- frequency analysis e c (u) if we write 

/•OC 

e c = I e c {u)d{2mu/hk 2 F ) . (23) 
Jo 

Gori-Giorgi et alO& have obtained an "exact" wave- vector analysis through the Fourier transform of the coupling- 
averaged correlation-hole density and built an analytic model for it by combining exactly known limiting behaviors 
with Diffusion Monte Carlo data, as in the work of Perdew and WangSl on the 3D EG. 

Figure 2 compares the results for e c (q) given by the theories of present interest with the QMC ones that have kindly 
been communicated to us by Dr. P. Gori-GiorgiS 6 .. While all theories perform similarly well in this test, it is pleasant 
to notice the reasonably good behavior shown by the ALDA. Finally, the imaginary-frequency analysis of the various 
theories is reported in Figures 3 and 4, in the low-u and large-it regimes respectively. As in the 3D case^, in all 
theories s c (u) starts with a finite negative value at it = and vanishes at large u. In ALDA e c {u) becomes positive 
at u ~ 7hk F /(2m) and ultimately vanishes from above. 

To conclude our work, in Appendix A we briefly comment on how the coupling-constant dependence of the frequency- 
dependent exchange-correlation kernel in an inhomogeneous electronic system of arbitrary dimensionality may be 
found from the knowledge of its density dependence at full coupling strength. 
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APPENDIX A: SCALING OF THE INHOMOGENEOUS EXCHANGE-CORRELATION KERNEL 

From time-dependent density functional theory and dynamical scaling Lein et al& have derived a coordinate-scaling 
relation in dimensionality D=3 for the coupling-constant dependence of the exchange-correlation kernel f xc [n] (r, r'; uS) 
of an inhomogeneous many-electron system. 

It is easy to generalize their proof to an inhomogeneous system in arbitrary space dimensionality D: namely, it is 
possible to prove that 

/ x A c M(r,r';c) = A 2 /xc[n'](Ar, Ar'; u;/X 2 ) (24) 

where n'(r, t) — A _D rt(r/A, t/X 2 ). In fact, the dimensionality of the system enters the proof only insofar as integrations 
of the time-dependent density n(r, t) over space are involved. 
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FIG. 1: Top panel: deviation of approximate correlation energy from the QMC correlation energy of Rapisarda and Senatore 
(in units of Ryd) as a function of r s for 0.5 < r s < 10 (the AKA results are omitted for r s 8 where they do not seem to vary 
smoothly). Bottom panel: approximate and QMC correlation energy (in units of Ryd) as a function of r s for < r s < 0.8 (the 
QMC results also embody the limiting form in Eq. (21)). 





FIG. 3: Imaginary-frequency analysis of the correlation energy per electron (in units of Ryd) at r 3 — 1 as a function of 
2mu/(hk%) in the low-u regime. In the inset an enlargment of the results from ALDA, QSTLS and DPGT is shown. 
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FIG. 4: Imaginary-frequency analysis of the correlation energy per electron (in units of Ryd) at r s 
2mu/(hkp) in the large-it regime. 
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